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^ ' Short-distance QCD is employed to calculate the J/ ip survival prob- 

ability in an equilibrating parton gas, whose evolution is governed by a 
^ ' set of master rate equations. Partons in the early stage of high-energy 

nuclear collisions may initially not be in equilibrium, but their average 
transverse momentum is sufficiently high to break up a QQ bound state. 
Such a break-up during the evolution of the parton gas is shown to cause 
a substantial J /if) suppression at both RHIC and LHC energies, using 
realistic estimates of the initial parton densities. The transverse mo- 
mentum dependence of the suppression is also shown to be sensitive to 
the initial conditions and the evolution history of the parton plasma. 
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I. INTRODUCTION 



It is generally believed, and confirmed by lattice QCD calculations 0, that 
hadronic matter under extreme conditions will form a plasma in which quarks and 
gluons are no longer confined to individual hadrons and are in both thermal and 
chemical equilibrium. To search for such a quark-gluon plasma, it is proposed to 
study collisions of heavy nuclei at extremely high energies. Recent work using models 
based on perturbative QCD indeed shows that at high energy a dense partonic system 
can be produced Though it is not all clear if such partonic systems will reach 

thermal and chemical equilibrium before hadronization |J7|— |TT[| , the partons, mainly 
gluons, are certainly in a deconfined state at such high densities |12[| . 

Many signals could arise from this deconfined state, such as charm quark enhance- 
ment fl3Hl"5f or enhanced photon and dilepton production |16,17|. Charm quarks, for 
example, cannot be easily produced during the mixed and hadronic phases of the 
dense matter, due to their large masses compared to the temperature. They can be 
readily produced only during the early stage of the evolution, when partonic degrees 
of freedom are relevant. In this paper, we will discuss how pre-equilibrium J/ip sup- 
pression can be used to probe the early deconfined state of the partonic system and 
the dynamics governing its evolution toward equilibrium. 

J/ip suppression due to color screening has been proposed to probe deconfine- 
ment [[OJ. This requires that the interactions of J/ip with hadrons and deconfined 
partons are different Because of its small size, a heavy quarkonium can probe the 
short-distance properties of light hadrons. It is thus possible to make a parton-based 
calculation of J/^-hadron cross section via an operator product expansion method 
similar to that used in deeply- inelastic lepton-hadron scatterings JT^|-|2l| . The result- 
ing J/^-hadron cross section can be related to the distribution function of gluons 
inside a hadron. The energy dependence of the cross section near the threshold of 
the break-up of a J/ip is determined by the large x behavior of the gluon distribution 
function, giving rise to a very small break-up cross section at low energies. Only at 
very high energies, this cross section will reach its asymptotic value of a few mb. In 
other words, the dissociation can only occur if the gluon from the light hadron wave 
function is hard in the J ftp's rest frame, i.e. its energy is high enough to overcome 
the binding energy threshold. A hadron gas with temperature below 0.5 GeV cer- 
tainly cannot provide such energetic gluons to break up the J/ip. Therefore, a slow 
J/ip is very unlikely to be absorbed inside a hadron gas of reasonable temperature 
|T9|. This conclusion does not seem to be affected substantially by non-perturbative 
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effects, analyzed in Ref. . 

On the other hand, a deconfined partonic system contains much harder gluons 
which can easily break np a, J/ip [pj3,E3|]. A study of the energy dependence of the 



gluon- J/ip inelastic cross section [53 shows a strong peak just above the break-up 
threshold of the gluon energy, eo = 2Mb — Mj/^, where Mj/^ and Mo are the 
J/ip and D meson masses, respectively In the pre-equilibrium stage, i.e., before 
the partons have reached equilibrium, the average parton transverse momentum is 
sufficiently large |24|] to break up a J/ip, provided the partons are deconfined. The 
dissociation of J/ip will continue during the whole equilibration process until the 
effective temperature drops below a certain value or the beginning of hadronization, 
whichever takes place first. Therefore measurements of J/ip suppression can probe 
the deconfinement of the early partonic system and shed light on the subsequent 
equilibration process, provided that possible nuclear effects on the production of QQ 
pairs and on pre- resonance charmonium states are understood and taken into account. 

In the following we will first calculate the thermal gluon- J/ip dissociation cross 
section at different temperatures and for different J/ip transverse momentum. We 
then follow the evolution of an initially produced parton gas toward equilibrium and 
calculate the resulting total survival probability of a J/ip and its P^-dependence. 



II. J/ip DISSOCIATION BY GLUONS 

The operator product expansion allows one to express the hadron-J/-?/> inelastic 
cross section in terms of the convolution of the gluon- J/ip dissociation cross section 



with the gluon distribution inside the hadron fL9 |. The gluon- J/ip dissociation cross 
section is given by [^3 



^2^/32^/^ 1 (g°Ao-l) 3/2 



3 V3; \3g*J ml {q^/e.f ' 

where g s is the coupling constant of gluon and c quark, uiq the c quark mass, and 
q° the gluon energy in the J/ip rest frame; its value must be larger than the J/ip 
binding energy eo- Since for the tightly bound ground state of quarkonium the binding 
force between the heavy quark and antiquark is well approximated by the one-gluon- 
exchange Coulomb potential, the QQ bound state is hydrogen-like and the Coulomb 
relation holds, 




e o = T^t m Q ■ (2) 
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The cross section thus can be rewritten as 



2 V 32 V 1 fa°Ao-D 3/2 (3) 

W ' 3 U/ rngfeomg)^ ( 9 o/ £o) 5 ■ I > 



As shown in Monte Carlo simulations ||24|| , the parton density in the early stage of 
high-energy heavy-ion collisions has an approximate Bjorken-type |25j scaling behav- 
ior. We will only consider J/i[) suppression in the central rapidity region (yjM — 0). 
In this case, the J/ip will move in the transverse direction with a four- velocity 



u 



(M T ,P T ,0)/Mj^ (4) 



where Mr = \jPt + ^-j/ip * s defined as the J/ip's transverse mass. A gluon with a 
four-momentum k = (k°, k) in the rest frame of the parton gas has an energy q° = k-u 
in the rest frame of the J/ip. The thermal gluon- J/ip dissociation cross section is then 
defined as 

J d 3 kv re ia(k ■ u)f(k°;T) 
■ «)>* = S<Pkf(k°;T) ' (5) 

where the gluon distribution in the rest frame of the parton gas is defined as 

f^T) = ^— x (6) 

with X g < 1 specifying the deviation of the system from chemical equilibrium. The 
relative velocity v Te \ between the J/ ip and a gluon is 

Vrcl k°M T k°M T ■ [ ' 

Changing the variable to the gluon momentum, q = (q°,q), in the rest frame of the 
J/ip, the integral in the numerator of Eq. ([5]) can be rewritten as 



/ 



<i 3 9^(9°)/(*°;r), 



where 

k = (q°M T + q-P T )/Mj^. (9) 

One can carry out the integral in the denominator, / d 3 kf(k°;T) = 8ir((3)X g T 3 , and 
the angular part in the numerator, to get 

M w , ,3/2^ r ,>-i) 3/2 , 



X 4 



3 C(3)irM T T d m Q ^ A 

(10) 
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with T n = T/n and 



± e M T ±P T 
a n = 7p T7 • (11) 



In order to understand the temperature and Pt dependence of the thermal gluon- 
J/if) dissociation cross section, we first plot in Fig. 1 the cross section a(q°) of Eq. (§) 
as a function of the gluon energy in the J/ip' rest frame. It decreases strongly toward 
the threshold and is broadly peaked around q° = 10eo/7 = 0.92 GeV, with a maximum 
value of about 3 mb. Low-momentum gluons have neither the resolution to distinguish 
the heavy constituent quarks nor the energy to excite them to the continuum. On 
the other end, high- momentum gluons also have small cross section with a J/ip since 
they cannot see the large size. 

We can also express the cross section as a function of the center-of-mass energy of 
gluons and the J '/if), cr(q°) = cr(s/2Mj/^ — Mj/^/2), where s = (k + Pj/^) 2 . One can 
thus translate the energy dependence in Fig. 1 into temperature and Pt dependences 
after thermal average, since the thermally averaged < s > is proportional to both Pt 
and temperature T. In Fig. 2 we plot the thermally averaged gluon- J/if> dissociation 
cross section as a function of temperature for different values of the J/^'s trans- 
verse momentum Pt- We observe the same kind of peak structure, with a decreased 
maximum value due to the thermal average. The position of the peak also shifts 
to smaller values of T when Pt is increased, corresponding to a fixed value of the 
averaged center-of-mass energy (s). A similar behavior is expected if one plots the 
thermal cross section as a function of Pt at different temperatures, as done in Fig. 3. 
However, in this case, the peak simply disappears at high enough temperatures, be- 
cause the averaged (s) will be above the threshold value even for Pt = 0. These 
features will have considerable consequences for the survival probability of a J/ip in 
an equilibrating parton gas, especially the Pt dependence. We should also mention 
that the use of the Bose-Einstein distribution function has an effect of about 20% on 
the thermal cross section, compared to that obtained with a Boltzmann distribution 

m 



III. J/ip SUPPRESSION IN AN EQUILIBRATING PARTON GAS 

Using the thermal cross section just obtained, we can now calculate the survival 
probability of J/ip in an equilibrating parton plasma. In this paper, we will neglect 
the transverse expansion and consider only longitudinal expansion. We will also only 
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consider J /if) suppression in the central rapidity region. A J /if) produced at point f 
with velocity v in the transverse direction will travel a distance 



d= -rcos(f) + ^R 2 A -r 2 (l -cos 2 0) (12) 

in the time interval = Mt(1/Pt before it escapes from a gluon gas of transverse 
extension Ra; here cos0 = v ■ f. Suppose the system evolves in a deconfmed state 
until the temperature drops below a certain value, which we assume to be 200 MeV. 
The total amount of time the J/if> remains inside a deconfmed parton gas is the 
smaller one of the two times t$ and tf, the life-time of the parton gas. Assume 
that the initial production rate of the J /if} is proportional to the number of binary 
nucleon-nucleon interactions at impact-parameter r, N A (r) = A 2 (l — r 2 j R a )/2tcR\. 
The survival probability of the J / if) averaged over its initial position and direction in 
an equilibrating parton gas is 

q/ u \ I d 2 r[R\ - r 2 ) exp [- drn g (r)(v rel a{k ■ u)) k ] 

S{Pt) ~ Jd 2 r(R 2 A -r 2 ) ' (13) 

where 

t min = min(t^,t f ) , (14) 

and n g (r) is the gluon number density at a given time r. 

In Eq. (|I~3D, both the gluon number density n g (r) and the thermal cross section 
(v re \a(k ■ u))k depend on the temperature, which in turn is a function of time. In ad- 
dition, n g is proportional to gluon fugacity which also evolves with time. To evaluate 
the survival probability we need to know the entire evolution history of the parton 
system. Roughly speaking, one can divide this history into two stages: (1) First there 
is kinetic thermalization, mainly through elastic scatterings and expansion. The kine- 
matic separation of free-streaming partons gives us tq ~ 0.5-0.7 fm/c as an estimate 
of the time when local isotropy in momentum distribution is reached PJ2"1|. (2) The 



parton gas now further evolves toward chemical equilibrium through parton prolifera- 
tion and gluon fusion; it does so until hadronization or freeze-out, whichever happens 
first. This evolution can be determined by a set of master rate equations which give 
us the time dependence of the temperature and fugacities. 

In this paper we will not address the question of pre- resonance J /if) suppression, 
which has been discussed and shown to be responsible for the J /if) suppression ob- 



served in p — A and S — U collisions |26|- p9fl . We will here consider the suppression of 
fully formed physical J /if) states as it should take place if nuclear collision produces 
dense partonic system. 
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Following Ref. ||, we characterize the non-equilibrium of the system by gluon 
and quark fugacities which are less than unity. The dominant reactions leading to 
chemical equilibrium are assumed to be the following two processes: 

99 «-> 999, gg^qq ■ (15) 

Assuming that elastic parton scatterings are sufficiently rapid to maintain local ther- 
mal equilibrium, the evolution of the parton densities can be given by the master 
rate equations. Combining these master equations together with one-dimensional 
hydrodynamic equation, one can get the following set of equations 0: 

^ + 3^ + - = Rs(l - X g ) - 2R 2 (1 - §), (16) 

£ + si + ± = (17) 

Ag 1 T bl Ag Ag 

(Ag + — A g ) 3/4 T 3 r = const., (18) 

&2 



where a x = 16C(3)/vr 2 « 1.95, a 2 = 8tt 2 /15 » 5.26, b x = %{3)N f /ix 2 « 2.20 and 



b 2 = 7ir 2 Nf/20 « 6.9. The density and velocity weighted reaction rates 



^3 = \{Ogg^gggV)n gi R 2 = \{^gg-,qqV)n g , (19) 



can be found in Refs. PJI5|. After taking into account parton screening jTJJ and 



Landau- Pomeranchuck-Migdal effect in induced gluon radiation [p0|,31]] , R3/T and 



R 2 /T are found to be functions of only A g . Solving the above rate equations as shown 
in Ref. 0, one finds that the parton gas cools considerably faster than predicted 
by Bjorken's scaling solution (T 3 r = const.), because the production of additional 
partons approaching the chemical equilibrium state consumes an appreciable amount 
of energy. The accelerated cooling, in turn, slows down the chemical equilibration 
process. For the initial conditions given by HIJING Monte Carlo simulations [Q, 
at RHIC energy the parton system can hardly reach its equilibrium state, since the 
effective temperature here drops below T c m 200 MeV in a time of some 1-2 fm/c. 
At LHC energy, however, the parton gas comes very close to equilibrium, since the 
system may exist in a deconfined state for as long as 4-5 fm/c. In the following we 
will use the numerical results for the time evolution of the temperature and fugacities 
obtained from the above master equations to calculate the survival probability of a 
J/ip in such an equilibrating parton system. 

Shown in Fig. 4(a) are the J/ip survival probabilities in the deconfined and equi- 
librating parton plasma at RHIC and LHC energies with initial conditions given by 
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HIJING Monte Carlo simulations ||J], denoted as set 1 in Table |. We find that there 
is stronger J/ip suppression at LHC than at RHIC energy, due both to the higher 
initial parton densities and longer life-time of the parton plasma. The increase of 
the survival probabilities with J/ip's transverse momentum is a consequence of the 
decrease of the thermal cross section with increasing Pt at high temperatures, as 
shown in Fig. 3, and the shorter time spent by a higher- Pj- J/ip inside the parton 
plasma, an effect first considered in Ref. P^] . J/ip dissociation during the late stage 
of the evolution should have a peak and increase a little with Pt when the temper- 
ature drops below 0.3 GeV, as illustrated in Fig. 3. This behavior flattens the total 
survival probability at small values of Pt when we integrate over the entire history of 
the evolution from high initial temperatures. For a parton system with a low initial 
temperature (below 300 MeV), the Pt dependence of the survival probability should 
be even flatter. One can therefore use the Pt dependence to shed light on the initial 
temperature and the evolution history of the system. 

To demonstrate the effects of the chemical non-equilibrium in the initial system, 
we also show in Fig. 4(b) the survival probabilities in an ideal parton gas with the 
same initial temperatures as before, but with full chemical equilibrium (unit initial 
fugacities). In this case, the temperature simply decreases like T(r) = T{tq){t g /t) 1 / z . 
We see that the J/ip is now much more suppressed than in the case of an equilibrating 
parton plasma, because of both the higher parton density and the longer life time of 
the system. 

Since there is considerable uncertainty in the estimate of the initial parton pro- 
duction by HIJING Monte Carlo simulations, as discussed in Ref. 0, we would like 
to test the sensitivity of J / ip suppression to this. We therefore multiply initial parton 
number densities by a factor of 4, thus increasing initial parton fugacities. We denote 
such initial conditions as set 2 in Table |I|. If the uncertainties in initial conditions are 
caused by the soft parton production from the color mean fields, the initial effective 
temperature will decrease. Therefore, we can alternatively increase the initial parton 
density by a factor of 4 and at the same time decrease T Q to 0.4 and 0.72 GeV at 
RHIC and LHC energies, respectively. This leads to higher initial fugacities, listed as 
set 3 in Table |I[ The corresponding survival probabilities calculated with these two 
sets of initial conditions are shown in Fig. 5. We can see that the J/ip suppression 
is much stronger if the initial parton densities are higher. Comparing the solid and 
dashed lines, however, shows that the J / ip suppression is less sensitive to the variation 
of the initial temperature and fugacities as far as the parton densities are fixed. 
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IV. CONCLUSIONS 



To summarize, we have used the cross section of J/ip dissociation by gluons to 
calculate the J / ip suppression in an equilibrating parton gas produced in high-energy 
nuclear collisions. The large average momentum in the hot gluon gas enables glu- 
ons to break up the J/ip, while hadron matter at reasonable temperature does not 
provide sufficiently hard gluons. We find a substantial J/ip suppression in such a 
non-equilibrium partonic medium; however, it is smaller than that in a fully equili- 
brated parton plasma. In particular, in an equilibrating plasma the behavior of the 
J/ip-gluon cross section at high gluon momenta reduces the J/ip suppression at large 
P T . 

In addition to the J / ip dissociation during the equilibration of the parton plasma, 
there are other possible sources of suppression for the actually observed J ftp's. As 
already noted, nuclear modifications of the QQ production process e.g., through mod- 



ified gluon distributions in a nucleus [[33|,|4[], multiple scattering accompanied by en- 



ergy loss |35|], or a suppression of the nascent J/ip before it forms an actual physical 
29fl must be taken into account. Such effects would cause J/ip suppres- 



resonance 



sion in addition to what we have obtained from the equilibrating parton plasma and 
modify the transverse momentum dependence of J/ip suppression [^6|. Moreover, 
interactions between gluons and cc bound states before the kinetic thermalization of 
the partons could also lead to a substantial J/ip suppression; this would depend on 
the time needed to achieve local momentum isotropy. On the other hand, gluon fusion 
could also result in J / ip production during the evolution of the parton system, similar 
to the enhancement of open charm. Although studies of pre-equilibrium open charm 
production indicate | JT5|JT5| | that J/ip production during the parton evolution is not 
significant compared to primary production, a consistent study of J/ip suppression 
should include the pre-equilibrium production in a form of a master rate equation. 
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TABLES 







LiiU^i ) 




LHC(2) 




LML-(3 j 


T(GeV) 


0.55 


0.82 


0.55 


0.82 


0.4 


0.72 




0.05 


0.124 


0.2 


0.496 


0.53 


0.761 


\ 


0.008 


0.02 


0.032 


0.08 


0.083 


0.118 


n 9 (fm~3) 


2.15 


18 


8.6 


72 


8.6 


72 


n g (fm _ 3) 


0.19 


1.573 


0.76 


6.29 


0.76 


6.29 



TABLE I. Different sets of initial conditions of the temperature, fugacities and parton 
number densities at tq = 0.7fm/c for RHIC and tq = 0.5fm/c for LHC. 
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FIGURES 

FIG. 1. Gluon- J/^ dissociation cross section as a function of the gluon energy q° in the 
rest frame of the J /ip. 

FIG. 2. The thermal-averaged gluon- J /ip dissociation cross section < v TC \a > as a func- 
tion of the temperature at different transverse momenta Pt ■ 

FIG. 3. The thermal-averaged gluon- J /ip dissociation cross section < v re \a > as a func- 
tion of the transverse momentum Pt at different temperatures. 

FIG. 4. (a) The survival probability of J/ip in an equilibrating parton plasma at RHIC 
and LHC energies with initial conditions given as set 1 in Table |, (b) and for an initially 
equilibrated plasma at the same temperatures. 

FIG. 5. The survival probability of J/ip in an equilibrating parton plasma with initial 
conditions given as set 2 (solid) and set 3 (dashed) in Table |. 
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